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Phases of a bilayer Fermi gas 
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We investigate a two-species Fermi gas in which one species is confined in two parallel layers and 
interacts with the other species in the three-dimensional space by a tunable short-range interaction. 
Based on the controlled weak coupling analysis and the exact three-body calculation, we show that 
the system has a rich phase diagram in the plane of the effective scattering length and the layer 
separation. Resulting phases include an interlayer s-wave pairing, an intralayer p-wave pairing, 
a dimer Bose-Einstein condensation, and a Fermi gas of stable Efimov-like trimers. Our system 
provides a widely applicable scheme to induce long-range interlayer correlations in ultracold atoms. 

PACS numbers: 03.75.Ss, 05.30.Flc, 74.20.Rp, 74.78.-w 



I. INTRODUCTION 

One of the central themes in ultracold atoms is to 
provide highly tunable model systems for other sub- 
fields in physics. A prominent example is the realiza- 
tion of ultracold two-dimensional atomic gases which 
have enabled detailed microscopic studies of Berezinskii- 
Kosterlitz-Thouless physics relevant to a wide variety of 
two-dimensional phenomena [1-3]. In addition to sim- 
ple two-dimensional systems with single layers, what 
have attracted considerable attention in condensed mat- 
ter physics are bilayer or multilayer systems. Here, extra 
degrees of freedom generated by layers and the long-range 
Coulomb interaction between them are expected to lead 
to intriguing physics such as an interlayer exciton con- 
densation in bilayer semiconductors [4, 5], quantum Hall 
bilayers [6], and bilayer graphenes [7]. 

In ultracold atoms, the analogous multilayer geometry 
can be created by confining atoms by a strong optical 
lattice in one direction. However, since the long-range 
Coulomb interaction is absent in neutral atoms, sepa- 
rated layers are simply decoupled without interlayer tun- 
neling. A novel scheme that we propose in this article 
to realize the long-range interlayer correlation in ultra- 
cold atoms is to use a mixture of two atomic species A 
and B and confine only A atoms in the multilayer ge- 
ometry with keeping B atoms in the three-dimensional 
space. Here the correlation between A atoms confined in 
different layers can be induced through the interaction 
with the background B atoms that are free to propagate 
from layer to layer. The advantage of this scheme is that 
one can tune the A-B interaction by interspecies Fesh- 
bach resonances and thus rich phenomena are expected 
to occur as a function of the interaction strength. Fur- 
thermore, such a system may be thought of as a tunable 
model system for a slab phase of nuclear matter by iden- 
tifying A (B) atoms with protons (neutrons) [8, 9]. 

In this article, we give detailed analyses in the case of 
a bilayer Fermi gas that can be realized by using a Fermi- 
Fermi mixture of, for example, ^Li and [10-12] and 
confining only one species in two parallel layers. Based 
on the controlled weak coupling analysis and the exact 
three-body calculation, we show that the system exhibits 



at least four distinct quantum phases as summarized in 
Fig. 1. In the weak coupling region, an effective attrac- 
tion between two A atoms mediated by B atoms leads to 
superffuidity due to an interlayer s-wave pairing (lower 
left phase) or an intralayer p-wave pairing (upper left 
phase) depending on the layer separation. On the other 
hand, in the strong coupling region, an A atom captures a 
B atom to form a tightly bound molecule and the ground 
state becomes a dimer Bose-Einstein condensation (right 
phase) . Finally, when the A-B interaction is close to the 
resonance, we find that two A atoms confined in different 
layers with one B atom form a three-body bound state 
leading to a Fermi gas of trimers in a dilute system (lower 
middle phase). Remarkably, there is an infinite number 
of such three-body bound levels exactly at the resonance 
resembling the Efimov effect in three dimensions [13]. 



The system under consideration is described by the 
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FIG. 1: Conjectured phase diagram at zero temperature in 
the plane of the effective scattering length (aoftfcp)"^ and the 
layer separation fcpd. Here, fep ~ fcpA ~ kvB is assumed. 
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FIG. 2: Effective interaction between two A atoms (solid line) 
induced by the interaction with B atoms (dotted line) in the 
weak coupling limit. Here, i,j — 1,2 are layer indices. 



action (here and below h = 1 and /cb = 1) 
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Here, tpAiitjX) with cc = (2;,?/) is a spinless fermionic 
field for A atoms confined in two parallel layers that are 
labeled by the index i = 1,2 and located at z — Zi. 
ipsityX, z) is another spinless fermionic field for B atoms 
in the three-dimensional space. The last term in the ac- 
tion describes the short-range A-B interaction, gg is a 
cutoff (A) dependent bare coupling, which can be elim- 
inated by introducing the physical parameter, the effec- 
tive scattering length Ooff, through ^ — V™b™-^b ^ _ 
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27ra„ff ■^^^^ = mArriB/imA + ms) being the 

reduced mass [14, 15]. Uctf is arbitrarily tunable by means 
of the interspecies Feshbach resonance and the limit 
tteff — > — (+)0 corresponds to the weak (strong) attrac- 
tion between A and B atoms. Throughout this article, 
the interlayer tunneling is assumed to be negligible and 
Fermi momenta of A and B atoms are defined through 
their densities; fcp^ = (47rn^)-^/^ and fcps = {Gn^nB)^^^- 



II. WEAK COUPLING LIMIT 

In order to elucidate phases appearing in our system, 
we start with the weak coupling limit Ucs — > —0 in which 
a controlled perturbative analysis is possible. Unlike the 
ordinary BCS-BEC crossover purely in two or three di- 
mensions, the A-B pairing at weak coupling does not take 
place in our mixed dimensional system because of the ab- 
sence of the Cooper instability [15]. What can happen 
instead are pairings between A atoms using the effective 
attraction induced by the interaction with the Fermi sea 
of B atoms. To the leading order in UeS, the back-to-back 
scattering of A atoms in layers i and j is described by 
the Feynman diagram in Fig. 2. The resulting induced 



interaction Vij{\p — q\) has rather a simple form in the 
real space: 
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where = y^x'^ + [zi — ZjY is the distance between 
the two A atoms. The oscillatory decaying factor is well 
known in the Ruderman-Kittel-Kasuya-Yosida interac- 
tion [16--18] and leads to the following physical interpre- 
tation of I^j (|a;|): The density modulation of background 
B atoms produced by one A atom in the layer i mediates 
the long-range interaction with the other A atom in the 
layer j. 

Now the pairings of A atoms are described by the BCS- 
type Hamiltonian: 
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In the mean-field approximation, the gap function is 
given by A,j{p) = nEq^^idP" Q\){i'AJ{-q)i'A^iq)) 
and the Fermi statistics of A atoms implies Aji{—p) = 
—A.ij{p). Although A atoms are spinless fermions, the 
layer indices i,j = 1,2 make more than one pairing 
pattern possible. When the layer separation is large 
d = |zi — Z2I — > 00, the interlayer interaction Vi^^j is 
suppressed while the intralayer interaction Vi^j is un- 
affected by d. Therefore, in this limit, intralayer pair- 
ings with Aii,A22 ^ are favored while their pairing 
symmetry has to be p wave [15]. On the other hand, 
when the layer separation becomes small d — >■ 0, Vi^j is 
no longer suppressed and thus an interlayer pairing with 
A12 = — A21 ^ is favored because the pairing symme- 
try can be s wave only in this singlet channel. Therefore, 
there has to be a quantum phase transition as a function 
of the layer separation. 

The quantum phase transition can be located by com- 
paring the energy densities for the interlayer s-wave 
and intralayer p-wave pairings. In the weak coupling 
limit, their energy densities are, respectively, given by 

(T^ind) = - ^^^^^ and (Hind) = - 

'"'^g'^"^ — ™d^2iL^ where the gap functions at |p| = kpA 
are given by Ai2/eFA oc e^'^/[™'*^i2'l and An/epA = 
A22/eR4 oc (p:^ ±ipy)e2'^/["-*^ii'l for negative V^f . 
Here, sfa = kpj^/{2mA) is the Fermi energy and v}^ = 

^ cos{l(p)Vij{\p — q\) with cosip = p ■ q is the par- 
tial wave projection of the induced interaction in which 
both incoming and outgoing momenta are restricted on 
the Fermi surface of A atoms; |p| — \q\ — kpA- It is 
the phase with the smaller energy density and hence the 
larger pairing gap that is realized in our system. 

The critical layer separation dividing the above two 
phases is plotted in Fig. 3 as a function of k-pA/k-pB, which 
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FIG. 4: Three-body scattering of two A atoms in different 
layers and one B atom. The double line represents the A-B 
scattering amplitude. 722 and 72i satisfy the same integral 
equations as 71 1 and 7i2. 



FIG. 3: (Color online) Critical layer separation fcpsd divid- 
ing the interlayer s-wave pairing (t^i2^ < t^i/') from the 
intralayer p-wave pairing (V^^^ > vii^) as a function of 

kpA/kFB- 



is independent of Ugs and mA/mB to the leading order in 

flcff. Below the critical separation where y/j^ < v}l\the 
interlayer s-wave pairing appears, while the intralayer p- 
wave pairing appears above the critical separation where 
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The two layers exhibit superfluidity in both 



phases and its critical temperature at weak coupling is 
the same order as the pairing gap at zero temperature; 
Tc ~ I Aij I . We note that our interlayer s-wave pairing 
bears an analogy to the interlayer exciton condensation 
in condensed matter systems [4-7]. Also it is worthwhile 
to point out that the intralayer px+ipy-wave pairing has 
a potential application to topological quantum computa- 
tion using vortices with non-Abelian statistics [19, 20]. 



tion in Eq. (1) has continuous symmetries of U(1)ai x 
U(l)yi2 X U(1)_B X 11(1)2 corresponding to particle num- 
ber conservations of A atoms in each layer, that of B 
atoms, and a rotation about z axis. An order parameter 
of the interlayer s-wave pairing {eijipAi4'Aj) 7^ breaks 
the full symmetries down to V{1)ai-A2 x U(1)b x U(1)z, 
while order parameters of the intralayer px+ipy-wave 

pairing {ipAiidx + idy)il)Ai) 7^ break them down to 
U(l)yii+yi2-z X U(1)b. Order parameters of the dimer 
Bose-Einstein condensation are {ip Aiip b) ^ {ip A21I' b) 7^ 0, 
which break the full symmetries down to \^{1)ai+A2-b x 
U(1)z. Therefore, the phases in the weak and strong cou- 
pling limits have to be divided by at least one quantum 
phase transition as a function of the interaction strength. 
Below we show that a novel phase can appear in between 
when the layer separation is small. 



III. STRONG COUPLING LIMIT 



IV. UNITARITY LIMIT 



We now turn to the strong coupling limit aos — +0, in 
which A atoms confined in layers capture B atoms from 
the bulk space to form two-body bound states (diniers) 
whose binding energy is given by E2 = —15 — ^-t-. The 

resulting system consists of the dimers localized around 
layers which interact weakly with each other and ex- 
cess B atoms. As long as the dimer size ~ aos is 
smaller than the layer separation d and the mean in- 
terparticle distance ~ fcp^, the dimers behave as two- 
dimensional bosons and therefore the ground state be- 
comes a dimer Bose-Einstein condensation in each layer. 
Accordingly, the two layers exhibit superfluidity up to 
the Berezinskii-Kosterlitz-Thouless temperature given by 



2-^11 



-In ^{-'^^'o-ndal^) in the limit a^s 



BKT ' 

[21], where M = niA + tub and Ud = ua are the 
dimer 's mass and density per layer. We note that the 
dimer Bose-Einstein condensates in different layers can 
interact with each other through the excess B atoms just 
as in Fig. 2 and Eq. (2). 

Although our system is found to exhibit superfluid- 
ity in both weak and strong coupling limits Ooff — >■ tOi 
their symmetry breaking patterns are different. The ac- 



The fourth phase realized in our system can be elu- 
cidated by studying a three-body problem of two A 
atoms conflned in different layers interacting with one 
B atom. Their scattering process is depicted in Fig. 4 
and all the relevant diagrams can be summed by solv- 
ing the integral equations for Tij{E;p,q). Here E is 
the total energy in the center-of-mass frame and p (q) 
is the relative momentum of the incoming (outgoing) A 
atom in the layer i (j) that does not scatter with the B 
atom flrst (last). The T-matrix elements have proper- 
ties 7ii — T22 and 7i2 = T21 and can be decomposed 
into even- and odd-parity parts under the exchange of 
one layer index; 7± = 7ii ± 7i2. After the partial wave 
projection T±{E]p,q) = ^ cos {I ip)T±iE;p,q) with 
cos (fi = p ■ q, we find T±^ to satisfy 



ri'^ (E; p, q) = TVmBTnAB if ^'^ (E + iO+;p, q) 
T2\e-p, k)K^^^E + iO+;k,q) 



dk k 
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FIG. 5: (Color online) Lowest binding energies of trimers for 
mass ratios mA/mB = 6.67 (bottom), 1 (middle), and 0.15 
(top) as a function of d/acB- The dashed and dotted lines are 
atom-dimer and dimer-dimer thresholds; E = E2 and 2E2. 



where K^^')(E;p,q) is given by 
I"' dip 



dJ -^^J p^+q^ + i^^pq cos ip-2mABE 



■ cos{lip) 



^pqcos(p-2mABE 



(5) 



The spectrum of three-body bound states (trimers) is 
obtained by poles of Ti''' as a function of E. When E ap- 
proaches one of the binding energies £'3 < — „ ^(°°»)^ ^ 

can write T^'^ asT2\E-p,q) ^ z'i\p,q)/{E - E^i). By 
solving the homogeneous integral equation from Eq. (4) 

satisfied by the residue ^ we find that the trimers 
exist only in an odd-parity (— ) and s-wave (/ = 0) chan- 
nel. The lowest binding energy e'^'^ as a function of 
d/ooff is plotted in Fig. 5 for three different values of 
the mass ratio m^/ms = 0.15, 1, and 6.67, correspond- 
ing to a mixture oi A — ^Li and B = ^°K, two dif- 
ferent internal states of the same atomic species, and 
A — '^^Y^ and B = ^Li. The trimer appears on the nega- 
tive side of the effective scattering length |aeff| ~ d and 
its binding energy decreases as d/a^a is increased. When 
d/cLcs +00, the trimer binding energy asymptotically 
approaches the atom-dimer threshold E = E2 from below 



as Es,/E2 - 1 cx e '^""VV -a »eff J ^ +0. 

Interestingly, when the A-B interaction is exactly at 
the two-body resonance |aoff| — > 00, there exists an 
infinite number of such trimer states whose spectrum 
is expressed by the form e'^^ — _ g-^Ti-Tt/so ^^^^^.^ for 
n — >■ 00 (not shown in Fig. 5). Such a geometric spec- 
trum at the resonance is well known as the Efimov effect 
in three dimensions [13]. Here the scaling exponent sq 
and the so-called Efimov parameter k can be determined 
as (so,k) = (0.741,0.0311), (0.828,0.0807), (1.30,0.231) 
for mA/mB = 0.15, 1, 6.67, respectively. An important 
characteristic of our trimer state in contrast to Efimov 
trimers in a free space is its stability against the three- 
body recombination because the two A atoms are spa- 
tially separated. 



The existence of stable AAB trimers, which has the 
same quantum number as eijipAiipAj'fpB, can lead to a 
Fermi gas of trimers in a dilute system. As long as 
the trimer size {2mAB\E3\)^^^^ is smaller than the 
mean interparticle distance ^ ^ea? the trimers behave 
as two-dimensional fermions localized around layers and 
form a Fermi gas with its density equal to ua- Although 
the trimer state exists even in the strong coupling limit 
Qoff +0, the trimer gas phase cannot persist there be- 
cause an addition of another B atom breaks up the trimer 
into two dimers when E3 > 2E2 is reached (see the dimer- 
dimer threshold in Fig. 5). Therefore, the trimer Fermi 
gas is realized only when the A-B interaction is close to 
the resonance jocff | > d. We note that AAB trimers with 
A atoms in the same layer and ABB trimers are absent 
for 0.0351 < mA/mB < 6.35 [14]. Once they are formed, 
the whole phase diagram would be dominated by such 
deeply bound trimers whose size is set by the thickness 
of layers. 



V. CONCLUSIONS 

We found that the dimer Bose-Einstein condensation 
appears in the strong coupling region where a^s ^ d and 
o-cskpA ^5 0(1), while the trimer Fermi gas appears in the 
unitarity region where |aoff| > d and kpAd ^ 0(1). As- 
suming that the rest of the phase diagram is occupied by 
the interlayer s-wave or intralayer p-wave superfluidity 
found in the weak coupling region, we conclude that the 
system has the rich phase diagram shown in Fig. 1. Our 
results can be tested by using a ultracold Fermi-Fermi 
mixture of ^Li and '"^K with a species-selective optical 
lattice. How the above four phases meet at the center of 
the phase diagram (shaded region in Fig. 1) and whether 
more quantum phases appear in our system are interest- 
ing open questions and can be addressed, in principle, by 
future experiments. In particular, further study on how 
the dilute gas of trimers evolves into a gas of their con- 
stituents as k-pAd is increased may shed light on a similar 
problem of the transition from nuclear matter to quark 
matter existing in the core of neutron stars. Finally, we 
emphasize that the scheme to induce long-range inter- 
layer correlations in ultracold atoms that is presented in 
this article by taking the bilayer Fermi gas as an example 
can be widely extended to multilayer geometries, mul- 
tiwire geometries, Bose-Bose mixtures, and Bose-Fermi 
mixtures. 
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